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ABSTRACT 

We present millimetre observations of a sample of 12 high redshift ultraluminous infrared 
galaxies (ULIRGs) in the Extended Growth Strip (EGS). These objects were initially selected 
on the basis of their observed mid-IR colours (0.0 < [3.6] - [4.5] < 0.4 and -0.7 < [3.6] - 
[8.0] < 0.5) to lie at high redshift 1.5 < z < 3, and subsequent 20-38/im mid-IR spectroscopy 
confirms that they lie in a narrow redshift window centered on z w 2. We detect 9/12 of the 
objects in our sample at high significance (> 3<j) with a mean 1200/im flux of < F^oopim > = 
1.6 ± 0.1 mJy. Our millimetre photometry, combined with existing far-IR photometry from 
the Far-IR Deep Extragalactic Legacy (FIDEL) Survey and accurate spectroscopic redshifts, 
places constraints both sides of the thermal dust peak. This allows us to estimate the dust 
properties, including the far-IR luminosity, dust temperature, and dust mass. We find that our 
sample is similar to other high-z and intermediate-z ULIRGs, and local systems, but has a 
different dust selection function than submillimeter-selected galaxies. Finally, we use existing 
20cm radio continuum imaging to test the far-IR/radio correlation at high redshift. We find 
that our sample is consistent with the local relation, implying little evolution. Furthermore, this 
suggests that our sample selection method is efficient at identifying ultraluminous, starburst- 
dominated systems within a very narrow redshift range centered at z ~ 2. 

Key words: galaxies: starburst - galaxies: formation - galaxies: high-redshift - submillimetre 
- radio continuum: galaxies - infrared: galaxies 



1 INTRODUCTION 



First discovered over 30 years ago, infrared (IR) luminous 
(LIRGs: L 8 -lOOOMm > 10 11 L Q ) and ultraluminous (ULIRGs: 
£8-iooo/Mn > 10 12 Lq) galaxies are among the most extreme ob- 
jects in the universe, with energy outp uts rivaling those of bright 
quasars (see Sanders & Mirabel 1996, for a r eview). Locally, the- 
oretical modeling (Mihos & Hernquist|[l99i) and both CO spec- 
troscopic (e.g.JSolomon et al.ll 19971) and optical/near-IR imaging 
dVeilleux et alj2002h indicate that they are powered in large part by 
star formation (SF) induced by major mergers. In the context of a 
merger-driven model of galaxy evolution, they represent the star- 
burst prelude to the rapid, self-regulated growth of a nuclear su- 
permassive black hole fe.gJSilk & Reesll 19981: [Murray et al]2005l ; 
|Pi Matteo et al.l2005l : lHopkins et alj2007l ; lYounger et alj2008ah . a 
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bright quasar, and eventually a "red and dead" early-type massive 
galaxy JSanders et al.lll988j|rllHopkins et al II 20061 . l2008bl fah ■ 

Though locally they contribute a very small fraction of 
the infrared luminosity density, at high redshift z > 1 LIRGs 
and ULIRGs take on increasing cosmol ogical importan c e and 
may dominate cosm ic SF at z > 2 telain et alj 1 19991 120021 : 
iLe Floc'h et al. I l2005h . There has also been recent observational 
evidence for a significant popul ations of hyperluminous sources 
at sti l l greater redshift z > 4 (Wang et al. 2007; Young er"et al.1 
l200ll2008bl : lDannerbauer et al.ll2008l ; lYounger et"alj|2008ch . Sev- 
eral techniques have been developed for selecting high redshift 
ULIRGs, includ ing: direct far-IR selection via (sub)millimeter 
surveys (e.g.. iHughes et all 1 1 9981: iBarger et ail 1 19981 : iGreve et ail 
120041 : IPope et al.l l2006l : IScott et all 120081) . dust-obscured galax- 
ies (DOGs) selected on the basis of their 24/im-R colour (e.g., 
iHouck et ai]|2005l ; lYan et al.l 12005 ), and those selected based on 
their mid-IR colours 1 Farrah et al] 20081 . Huang et al., in prepara- 
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Figure 1. The correlation between the millimetre-to-24^tm flux density ra- 
tio (_R(1200, 24 ) = Si200fj,m/S24fi) and the obscuration parameter of 
iLutz et alj i200^) : i_(24, 0.6) = logi (vF^(24/.m)/i;i^(0.6/.m)). Up- 
per limits are at the 3cr level are indicated by inverted, filled triangl es. We 
find the all of our sources fall above the mean measured bv lLutz et alj J2005L 
shaded region and dashed line), with a mean of < _R(1200, 24) >= 2.22± 
0.08 (dotted Une) with an intrinsic dispersion of c^( 12 oo 24) = 0.83. 
Therefore, our sources are more far-IR luminous at fixed 24fim indepen- 
dent of obscuration, which is consistent with a larger contribution from a 
starburst starburst over a dusty AGN. 



tion). Understanding the nature of these different populations, and 
the engine driving their extreme luminosities, is crucial to a thor- 
ough understanding of SF and galaxy evolution at high redshift. 

At the same time, a remarkably tight correlation between 
the radio and far-IR luminosity of galaxies, spanning several 
decad es in FIR luminosity, has been known for more than 20 
years faelou et alJl985l : ICondon & Broderickll 99 lMCondonl 19921 : 
Yun et al.ll20oTR 7 Though it is thought that this relationship arises 
from synchrotron loses associated with massive star for ming re- 
gions teurrimel|[l98r3 ; lHummel et al.ll 19881 ; [CondorJI 19921) . the de- 
tailed p hysics are relatively p oorly understood (see also the discus- 
sion in iThompson et al1l2006t) . Therefore, the persistence of this 
relation at intermediate to high redshift - where the conditions of 
galaxy formation and star forming environment were likely quite 
different than they are local l y - is a subject of grea t interest in cur- 
rent research ([Ga rrett 2002; iGruppioni et alj|2003 ; Appleton et al.1 
20041; iKoyacseTa l. 2006; Boyle et al.ll2007l : Iviahakis et alJl2007t 
Ibar et alfeOOlllSaiina et alj|2008l) . 

In this work, we present far-IR observations of a sample of 
12 mid-IR selected ULIRGs at z ~ 2, first identified by Huang 
et al. (in preparation). Our far-IR photometry, which spans both 
sides of the thermal peak, allows us to more robustly estimate the 
luminosity and dust properties - e.g., dust temperature and mass 
- of these objects. The relative uniformity of this sample, com- 
bined with mid-IR spectroscopic redshifts that indicate a remark- 
ably tight redshift range, makes it ideal for such a study. It further- 
more provides a promising platform for probing the far-IR/radio 
correlation at z ~ 2. 

This paper is organized as follows: in §[2] we outline our sam- 
ple selection, in §[3]we detail our observations, in §|4]we present the 
results of our observations, in §[5] we describe our far-IR spectral 
energy distribution fitting method, in §[6] we present estimates of 
the dust properties for objects in our sample, in §|7]we examine the 
far-IR/radio correlation at high redshift, and in §[8] we conclude. 
Throughout this work, we assume a flat concordance cosmology 



with f_ m = 0.3, Qa = 0-7 and h ■■ 
iTegmark et al.ll2004Ll2006l) . 



2 SAMPLE SELECTION 



0.7 dSpergel et al.l2003ll2007l; 



It has been shown dSawickil |2002| ; IHuang et all |2004| ; IPapovichI 
2008) that an efficient method for selecting high redshift galaxies is 
to track the 1.6/im stellar emission bump via mid-IR colors. Using 
the Infrared Array Camer a (IRAC: Fazio et al. 2004) on board the 
Spitzer Space Telescope ( I Werner et al.1 20041) . applying the color 
selection criteria 



0.0 < [3.6] - [4.5] < 0.4 and 
-0.7 < [3.6] - [8.0] < 0.5, 



(1) 

(2) 



where [3.6] denotes the AB magnitude in the 3.6/im band (and like- 
wise for the 4.5 and 8.0/im bands), will identify massive galax- 
ies in the redshift range 1.5 < z < 3.0 - the redshift range 
where the peak of the 1.6pm stellar bump is redshifted into the 
4.5 and 5.8^tm bands. Furthermore, since this is effectively a rest 
frame near-IR selection criterion, it is relatively insensitive to dust 
reddening and is roughly ste llar-mass limited (.Hua nget al.ll2004l 
120051; IConselice et alj2007ah . 

Huang et al. (in preparation) used this technique, combined 
with a 24/j.m flux limit of i^pm *5 0. 5 mjy usin g the Multiband 



Imaging Photometer for Spitzer (MIPS: lRieke et alj2004T> . to select 
high redshift, massive, IR luminous systems for mid -IR spectro- 
scopi c followup using the Infrared Spectrograph (IRS: Hou ck et al.1 
l2004h - both also on board Spitzer. This yielded a sample of 12 ob- 
jects in the Extended Groth Strip (EGS), all with Ltr ^ lO 12 !/© 
in a remarkably tight redshift range centered around z ~ 1.9. This 
was also a considerably more uniform sample of high redshift IR 
luminous galaxies than other Spitzer selected samples; see Huang 
et al. (in preparation) for a more detailed comparison to the var- 
ious other selection techniques. Our sample was drawn from the 
Huang et al. (in preparation) EGS sample with mid-IR derived 
spectroscopy redshifts. 



3 OBSERVATIONS 

The optical and near- to mid-IR obs e rvatio ns used for this work 
are described in detail in iDavis et al] | |2007|) . Huang et al. (2008, 
in preparation), and references therein. Briefly, we make extensive 
use of deep r' imaging with a limiting magnitude of 26.5 AB, IRAC 
imaging at (3.6, 4. 5, 5.8, 8.0)/xm to 5a flux limits of (1,1, 1.2, 6.3, 
6.9)pJy (see also iBarmbv et ai]|2004 l2008h. and MIP S imaging 
at 24/im to a 5<r depth of 77/iJy dSanders et alj 120071) . We also 
make use of extremely deep far-IR MIPS imaging at 70 and 160/mi 
that combined existing observations (PID 00008, PI Fazio) and the 
Far-Infrared Deep Extragalactic Survey (FIDELS; PID 30948, PI 
Dickinsonfl and deep 1.4 GHz radio imaging to a depth of 50/xJy 
dlvison et alj2007l) . Photometry for our sources was first presented 
by Huang et al. (2008, in preparation). 

The millimetre observations presented here were performed 
at the Institut de Radioastronomie Millimetrique (IRAM) 30 me- 
tre telescope at Pico Veleta during the Winter/Spring of 2008, us- 
ing the 117 element Max-Plank Bolometer Array 2 (MAMB02: 



1 For a description of the FIDELS observing strategy and data reduction, 
see http://data.spitzer.caltech.edu/popular/fidel/2007_sepl7/fidel_dr2.html 
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Table 1. The observed r' band fluxes and far-IR SEDs of the ULIRGs in our sample. Upper limits are at the 3<r level. The 70/im,160/im, and 20 cm photometry 
is from Huang et al. (in preparation). The 1200/tm MAMBO observations and data reduction are described in §|5] 




Figure 2. The millimetre detectability as a function of the 70/im- (i^Yo^m/ F\200^m'- left) ar| d 160/^m-to-millimetre (-Fi60Mm/-Fi200Mm: right) flux density 
ratio. Upper and lower limits are indicated by arrows. We find a strong negative correlation at the ~ 2<r level (see Table [2) between the millimetre flux and 
^70Mm/^i200Mm, reflecting the fact that warmer sources are less likely to have a significant cold dust component. We do not, however, find a significant 
correlation between the millimetre flux and -Fi60fim/^l200nm. 



iKrevsa et al.lll998t) . All our observations were performed in good 
weather conditions, with T230GHZ S 0.3 and 1.9 Hz (corresponding 
to a 0.52s wobbler "period") chopped sky noise < 100 mjy/beam 
(HPBW=11 arcsec). Primary pointing, focus, and flux calibration 
was performed on either Mars (~ 130 Jy) or another bright cali- 
brator source - including K3-50A (~ 7 Jy) and CW-LEO (~ 1.3 
Jy). Secondary pointing was performed on J1419+544 (~ 1 Jy), 
a bright quasar only 1 degree away. On-off observations were ob- 
tained in 20 minute scans with a wobbler throw of 35 arcseconds, 
which were repeated until either a significant detection (S/N > 3cr) 
or a noise level of ~ 0.4 mJy. This typically necessitated ~ 2 hours 
of integration time per source (full range 1-3 hours). The data were 
reduced using the MOPSIC Package developed by R. Zylka, using 
the standard reduction pipeline. 



4 RESULTS 

In all, 9 of the 12 targets had significant detections (> 3cr) at 
1200^im- egsl 1 was detected at ~ 2a, egs24 was detected at ~ la, 
and egsb2 was consistent with the sky noise. Excluding egsb2, the 
mean observed flux density < Fi 200pm >= 1.6 ± 0.1 mJy and 
mid- to far-IR colour < £(1200, 24) >=< S , i2oo M m/5'24 M m >= 



2.22 ± 0.08 with an intrinsic dispersion of a^ 



B(1200,24) 



0.83. 



The average 1200/im flux density and detection rate is similar to 
the IRAC selected starburst sample lLonsdale et al.l d2008l) . This av- 
erage flux density ratio noticeably lower than is typical of bright 
SMG s « £(1200,24) >^ 12: lEgami et ai]|2004l ; Hvison et al.1 
2004). However, while it is consistent with the detected sources 
in their sample, this is a significantly higher average value than 
has been observ ed in other Spit zer selected samples of high red- 
shift ULIRGs: dLutz et al.ll2005l) followed up a sample of 40 ob- 
scured 24^m sources (see also lHouck etai]|2005l : lYan et alj|2005l: 
IWeedman etal]|2006dlbl) with MAMBO at 1.2mm and found an 
overall < £(1200, 24) >= 0.74 ± 0.09 with a k(1200 2A) = 0.59, 
which is consistent with t he expectation from an IR bright ac- 
tive galactic nucleus (AGN: lElvis et al.lll994h . ISaiina et al.U2008h 
found similar millimetre fluxes for a comparable sample dYan etldl 
120071) . This is not unexpected, as obsc ured 24/xm sources ar e 
thought to contain a large fraction of AGN dWeedman et~afl2 006b ) . 

In Figure [T] we sh ow the obscuration pa- 
rameter from iLutz et al.l d2005l) . i?(24, 0.6) = 

logio {vF v (24/xm) /uF u (0.6/x m)), whe re vF v (0.6^tm) is 

taken from the r' imaging dDavis et al.l |2007|) . These authors 
found that when they restricted their sample to the most obscured 
sources with i?(24, 0.6) > 1.5, the average _R(1200, 24) increased 
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p 


Significance 


-FY0/jm/-Fl200Mm 


-0.76 


2.0 


-Fl60Aim / fl200/«n 


-0.10 


0.3 


-^7.7/jm 


-0.04 


0.1 
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0.26 


0.8 


EW 7 .7 M m 


-0.62 


1.6 


EWn.3 M m 


0.31 


0.7 



Table 2. The results of a Spearman rank correlation coefficient (p) analysis 
of the trends of different PAH parameters (see Figure [3] with the observed 
millimeter flux density -Fi200nm ■ 



to 0.74 ± 0.09 with cr^ (1200 24) = 0.59. Though they could 
not examine it directly due to the low detection rate of sources 
in their sample, this suggests a correlation between -R(24, 0.6) 
and i?(1200, 24). For our sources with S/N ^ 2, a Spearman 
rank correlation coefficient analysis finds a mild correlation with 
p = 0.27 at the ~ lcr level. The low significance of this correlation 
may be due to smal l sample size, but it is broadly consistent 
with the findings of iLutz et al. I d2005h and suggests that the most 
heavily obscured sources efficiently reprocess the attenuated stellar 
radiation into cold (< 60 K) thermal dust emission. 

Given the narrow redshift range of our sources, we can also 
use the 10pm and 160/im-to-millimetre flux density ratios to in- 
vestigate the correlation between strong millimetre emission and 
the dust properties of objects in our sample. At z ~ 2, the 70 and 
160/^m MIPS channels probe the rest-frame 24 and 50/^m emis- 
sion. These wavelengths probably probe AGN-heated dust (e.g., 
lElvis et alj [l994l; lUrrv & Padovanj [l995h and cold dust from a 
starburst ^Sanders & Mirabell 1996 ) respectively. In Figu 



s|2h 



show the millimetre flux as a function of both i 7 Vo Mm /^ ? i200/jm 
and -Fi60Mm/Ji200nm respectively, and in Table [2] we present 
the results of a Pearson rank correlation coefficient analysis. We 
find a strong negative correlation between the millimetre flux and 
-FYoMm/-Fi200(im; sources with more hot relative to cold dust - and 
therefore likely a larger AGN contribution - are less likely to be 
detected at 1200/itn. We do not however, find a strong correlation 
between Fieojun/.Fi200jun and Fi2oo M m- We believe that this is 
mostly likely a result of the relatively poor sensitivity of the \6Qpm 
FIDELS survey: if an object has enough cold dust to be detected at 
160/im at all, it is also likely to be a strong millimetre source. 



far-IR (Rigopoulou et al.lll999: 


Wu et al.ll2005l:lBrandl et all2006fc 


Calzetti et al. 2007; Desai et al. 


2007; Shi et al.l 20071), and there- 



fore we would expect strong PAH sources to show strong millime- 
tre emission. To investigate this for our sample, which contains all 
strong PAH emitters, in Figure[3]we show the correlations between 
the millimetre flux {F^oo^m) and PAH luminosities (Lr.r^.m and 
in. »nm) and equivalent widths (EW7.7,im and EWn.3 Mm ) at 7.7 
and 1 1.3/itn. The results of a Spearman rank correlation coefficient 
analysis are presented in Table|2] All of these results are low signif- 
icance, likely due to the effects of a relatively small and homoge- 
nous sample. They are broadly consistent - considering the low sig- 
nificance of the 1/7. 7 Mm inverse correlation - with observations of 
local systems that found a correlation between the PAH luminos- 
ity and ongoing star form ation (e.g.. IWu et al.ll20o3 : iBrandl et al.l 
2006; ICalzetti et al.ll2007t) . In particular, we find a positive corre- 
lation between 11.3/^m PAH luminosity and equivalent width and 
the millimetre flux at a 0.8 and 0.7<r level of significance respec- 



tively. Since the 11.3/i m feature is more sensitive to obscuration 
(Rigopoulo u et alj|l999l) , sources with strong observed features at 
this wavelength are likely to be more intense starbursts, and thus 
more likely to be detected in the far-IR. However, we find a sig- 
nificant negative correlation between the 1.1pm equivalent width 
an d J ? i2oonm at 1.6c r significance. This differs from the results 
of ISaiina et al. (2008) who find a positive correlation (p — 0.56) 
between F^oo^m and EW7.7 Mm at ~ 3<r confidence for a sam- 
ple o f high redshift obscured Spitzer selected ULIRGs dYan et all 
120071) . While the correlation seen for our sources is, again, low- 
significance, it suggests that the relationship between strong PAH 
emission and far-IR luminosity is somewhat different in their sam- 
ple - a reasonable result consid ering the far diffe rent mid- to far- 
IR properties of similar objects dLutz et al]|2005l) are significantly 
different from those presented here. 



5 SPECTRAL ENERGY DISTRIBUTION FITTING 

Our millimetre observations constrain the shape of the far-IR SED 
of our sources with photometry on both sides of the thermal dust 
peak. Given a choice of a simple far-IR SED parameterization, this 
allows us to characterize the dust properties in a sample of very 
luminous, high redshift systems from the far-IR directly, without 
invoking additional assumpti ons - for example ass uming the local 
far-IR/radio correlation (e.g., Chapm an et alj200^) . For this work, 
we use a single tempe rature greybody fitting form first introduced 
by iHildebrand (1983), in which the thermal dust spectrum is ap- 
proximated as F v = Q v B v {Td), where B v is the Planck function, 
Qv = Qo{v/vq)P is the dust emissivity with 1 < j3 < 2, and Td is 
the effective - or emission weighted - dust temperature. This yields 
a dust spectrum that for hv / kTd ^ 1 takes the form: 



,3+P 



F v oc 



exp(hv/kTd) 



(3) 



This is then typically paired with a power-law component at short 
wavelengths F v ~ v~ a , to account for a subdominant warmer dust 
component arising either from the starburst or an AGN, which is 
matched to the long wavelen gth SED at dF v /d v = — a to ensure 
a continuous first derivative. Blain et aT] d2003h find that this form 
provides a good match to the SEDs of observed systems from A ~ 
30 — 1000/im over a large range in total IR luminosities (see also 
Blainll 19991: ICarilri & YurJ[l 999; Dun ne et alj|2000l : lYun & Carillil 



2002l;lsaiina et alj2006l : lYang et alj2007 . for additional discussion 
of far-IR fitting methods) 

This model was fit to the 70, 160, and 1200/im data using 
a Levenberg-Marquardt least-squares minimization routine, the re- 
sults of which are summa rized in Figu res[4]and[5]and Table[3] Since 
T d is degenerate with /3 telain etalj r2003; Saii na et al.ll2006|) for 
sparsely sampled far-IR SEDs, we fix j3 — 1.5 which is consistent 
with the results of SED fitting of loca l and high redshift sy stems 
jDunne et alj2000l : lYun & Carillil2002l ; [Yang & PhiHipsl2007t) : us- 
ing a high er value of 13 will l ead to systematica lly lower dust tem- 
peratures jBlain et all 120031 : ISaiina et alj|2006t) . Furthermore, the 
accurate redshifts supplied by mid-IR spectroscopy were essen- 
tial in removing the additional degeneracy between Td and r edshift 
dCarilli & Yunfl 9991 : lYun & Carillill2002l ; iBlain et al.ll2003l) . From 
these fitted models we can estimate Td and the far-IR luminosity. 
We can also estimate the total dust mass in these systems according 
to: 



M d 



rest (Arest j Td) 



(4) 
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Figure 3. The correlation between the observed millimeter flux density i*i200^m and PAH luminosities (left: L7.7 Mm and Lii^m) and equivalent widths 
(right: EW7.7^ m and EWii.3 Mm ) at 7.7 (top) and 1 1.3/im (bottom). The results of a Pearson rank correlation coefficient analysis are shown in Table|5] 



where Md is the total dust mass, S u is the observed flux density, 
Dl is the luminosity distance, K(A res t) is the rest frame dust mass 
absorption coefficient at the ob served wavelength (for whi ch we 
use the Milky Way dust model o f Weingartn er"& Drainel2o6"lh . and 
B v (A ros t , Td) is the Planck function at the rest wavelength. At long 
wavelengths - e.g., 1200/im observed, ~ 400^im in the rest frame - 
observations are directly probing populations of dust that are pref- 
erentially colder than Td, which is effectively an emission weighted 
average over dust populations at a large range of physical tempera- 
tures. Therefore, we use the fitted models to estimate the peak flux 
and wavelength, which is dominated by dust at Td and simultane- 
ously minimizes the effects of uncertainty in (3 (k v oc v -13 ), for 
input into equation [4] 

Finally, we com pute the q paramete r for the far-IR/radio cor- 
relation as defined bv lHelou et al.l dl985l) : 



<1 = 



log ( 



i40-120/jm 

3.75 x 10 12 W 



log( 



il.4GHz 

W Hz- 



(5) 



where Z/4o-i20nm is the integrated IR luminosity between 40 and 
120/im, and Z/i.4gh z is the rest-frame 1.4 GHz continuum radio 
luminosity. Since our sources are at z ~ 2, their observed 1.4 GHz 
flux corresponds to a rest-frame frequency of 4 .2 GHz. Therefore , 
we applied a fc-correction similar to that used in lYang et al.l d2007l) . 
assuming a power-law synchrotron F„ jr adio oc v~ a " r , whe re we as- 
sume a r = 0.8 as in typical local star-forming systems dCondonl 
Il983lll992h . The results, as with the far-IR luminosity and dust 
properties, are listed in Table(3] and discussed in the following sec- 
tions. 



6 DUST PROPERTIES 

Our far-IR and millimetre photometry that constrains both sides of 
the cold dust peak allows a determination of the far-IR luminosity, 
dust temperature, and dust mass from the far-IR alone. As a result, 
the dust properties in our high redshift ULIRG sample can be deter- 
mined with a minimum of additional assumptions. While the prop- 
erties of dust in these extreme environments and at high redshift is 
of itself an interesting result, it is particularly so in the context of 
analogous and related objects at different redshifts and with differ- 
ent selection techniques. Therefore, we assembled results from five 
samples with extensive far-IR observations from the literature for 
comparison: another Spitzer-selected high redshift sample (HZS), 
a low redshift sample (LZS), intermediate redshift sample (IZS), 
a submillimetre selected sample (SMS), and finally a sample of 
quasars - which are thought to follow the UL IRG phase in an evo- 
lutiona ry sequence ( Sande rs et alJll988dlbTlct iHopkins et al]|2006l. 
2008b|a|) - with extensive far-IR observations (QS). Our far-IR 
SED parameterization is consistent with those used for the analysis 
of all samples discussed below - including the HZS, IZS, SMS, and 
QS - and therefore there are no hidden systematic biases. 

The HZS consists of bright, ob scured 24/im sources sected 
on th e basis of their -R(24, 6) ratio dHouck et al.ll2005l;lYan et al.l 
Weedman et al. 2006a|l bJ), wh ich have mid-IR spectrscopic 
dYan et alj|2007l : ISaiina et al.ll2007l) a nd MAMBO 1200 /xm pho- 
tometric foll owup dLutz et al.l l20ol : ISaiina et alj 120081) . As we 
note in § [4] iLutz et al] d2005l) find a signficantly lower average 
_R(1200, 24) ratio than the objects in our sample, even among the 
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Figure 4. The observed far-IR SEDs (see Table [T} of the z ~ 2 ULIRGs in our sample, along with the fitted model. The model parameters and derived 
quantities are listed in each panel, as well as tabulated in Table [3] and the fitting procedure summarized in §|5] The solid lines have a fixed dust emissivity 
j3 = 1.5, while the dotted lines are the same for = 2.0. 
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z 


T d 


i40-1000fim 


-£<40-120/jm 


M d 


9 






[K] 


[10 12 L o ] 


[1O 12 L ] 


[1O 8 M ] 




egsl 


1.88 


36 ±6 


2.1 ± 0.9 


1.7 ±0.7 


3.6 ± 1.6 


2.15 ±0.20 


egs4 


1.90 


37 ±3 


3.3 ±0.7 


2.7 ±0.6 


4.4 ±0.9 


2.29 ±0.11 


egslO 


1.94 


47 ±3 


5.5 ± 0.7 


5.0 ±0.7 


2.1 ±0.3 


2.44 ±0.08 


egsl 1 


1.80 


41 ±6 


2.5 ± 0.8 


2.2 ±0.7 


1.8 ±0.6 


2.25 ±0.15 


egsl 2 


2.01 


39 ± 3 


3.6 ±0.8 


3.1 ±0.8 


3.7 ±0.8 


2.57 ±0.15 


egsl4 


1.86 


41 ±1 


9.7 ±0.7 


8.3 ±0.6 


8.3 ±0.6 


2.21 ±0.04 


egs21 


3.01 


62 ± 5 


8.9 ± 1.5 


8.4 ±1.4 


0.8 ±0.1 


2.32 ±0.09 


egs23 


1.78 


44 ±3 


6.4 ±0.6 


5.6 ±0.6 


3.4 ±0.3 


2.43 ± 0.06 


egs24 


1.85 


43 ± 14 


1.2 ± 0.8 


1.1 ±0.7 


0.7 ±0.5 


2.07 ±0.30 


egs26 


1.80 


39 ±4 


2.6 ± 0.6 


2.2 ±0.5 


2.8 ±0.6 


2.10 ±0.11 


egs24a 


2.12 


34 ±4 


4.0 ± 1.5 


3.1 ± 1.1 


9.5 ±3.3 


2.03 ±0.15 


egsb2 


1.59 


34 ±3 


2.3 ± 0.6 


1.8 ±0.4 


5.6 ±1.3 


2.00 ±0.11 



Table 3. Derived parameters for our high redshift ULIRG sample, including: the mid-IR spectroscopic redshift (from Huang et al. in preparation) and the 
effective dust temperat ure (T^), far- IR luminosity from 40-1000^tm (-L40-iooonm ) and 40-120/xm (L40— i20um), dust mass (Mj), and far-IR/radio q 
parameter (see §|7]and Condon 1992). For details on the fitting procedure, see §|5]and Blain et aT] fe003l) . 



most obscured (-R(24, 0.6) > 1.5) subsample. The authors find a 
median dust temperature of 32 ± 8 K (35 ± 8 K for 14 strong-PAH 
srouces, and 32 ± 8 K for 34 weak-PAH sources), as compared 
to 41 ± 5 K for our sample of 12 objects, with comparable far- 
IR luminosities and dust masses. This difference may arise from 
the nature of their sources, combined with their SED fitting tech- 
nique; the HZS has sho wn a significant contribution to the far-IR 
by an AGN (see also S lTlWeedman et alj2006bl : ISaiina et al.l2008h . 



which the author s account for by fi tting a combination of models to 
the mid-IR (see lSaiina et alj|2007l) and far-IR separately in order 
to separate out the cold dust population. However, it is most likely 
driven largely by their choice of ft = 2, which leads t o cooler fit- 
ted dust temperatures th an the /3 = 1.5 adopted here dBlain et al] 
l2003l ; ISaiina et al.l2006t) - indeed, when we adopt (3 — 2.0 we find 
a median dust temperature of 33 ± 4 K while the total IR luminos- 
ity decreases by ^ 15%. Furthermore, since their mid-IR models 
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Figure 5. Same as Figure[4] 



will implicitly include some of the warm dust heated by the star- 
burst, their two-component fitting method may also lead to some- 
what lower dust temperatures. Therefore, while theirs is a rather 
different population of objects, we find that the dust properties in 
both are similar. This is an interesting resu lt, considering the large 
fraction of obscur ed AGN in such samples (Wee dman et alj2 006b; 
Sajina et al. 2008); the cold dust appears to be similar in both AGN 
and starburst dominated systems. 



The LZS is the SCU BA Local Universe Galaxy Survey 
(SLUGS lPunne et alj200(j|) . which consists of 104 local (z & 0.05) 
infrared luminous galaxies selected from the 60/xm flux limited 
IRAS Bright Source Catalog (BSC: ISoifer et alj[l98^) with com- 
plete 450 and 850^im followup with the Submillime tre Common 
User Bolometer Array (SCUBA: [Holland et al.l[l999l) at the J ames 
Clerk Maxwell Telescope (JCMT). T he IZS dYang et alj|2007l) was 
drawn from the FSC-FIRST Survey dStanford et alj|200d) . which 
consis ts of a matched sam ple of 60/im IRAS Faint Source Catalog 
(FSC: lMoshir et alj|l992h and 1.4 GHz sources from the Faint Im- 
ages o f the Radio Sky at Twenty Centimeters (FIRST: Bec ker et al.l 



1 19951) Sur vey, with complete 3 50^m followup from the SHARC— 
U camera dDowell et alj|2003l) at the Caltech Submillimetre Ob- 
servatory (CSO). The SMS dKovacs et alfeOOfj) consists of SMGs 
identified at 850/iin in blank field SCUBA surveys wit h 20cm ra- 
dio counterparts, optical redshifts dChapman et al J 20051) . and com- 
plete 350^tm follo wup with SHARCH -II (for a similar sampl e, see 
ICoppin et al]|2008l) . Finally, the QS dHaas et alj[2000l l2003h con- 
sists of 64 quasars sel ected from the Palomar Green Survey (PG: 
ISchmidt&Greenl 19831) with far-IR coverage from 5-200/^m from 
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Figure 6. The temperature-luminosity relation Lpjn — Tj, where Lpm 
is taken from 40-1000/im . The fitting method is detailed in §fj] and the 
values for our sample tabulated in Table fj] Incl uded are results fr om our 
sample (red circles) as compared to l ocal (D00: iDunne et alj|2000t black 
triangles), intermediate redshift (Y07: lYane et al.ll2007t open stars), high 



redshift submillimeter galaxies (K06: 



Kovacs et alj|2006l black squares), 



and PG quasars (H03 lHaasetai]|2003l blue circles). The dashed lines in- 
dicated tracks in constant dust mass (see Equation [7] and §[6) with (left to 
right) Md = 10 7 , 10 8 , 10 9 , and 10 10 Mq, assuming a fixed emissivity of 
j3 = 1.5, The dotted lines are the same for (3 = 2.0. 



ISOPHOT dLemkeet all 19961) on board the ISO satellite and some 
(sub)millimetre coverage with MAMBO and SCUBA. 

For single temperature dust models, the total far-IR luminos- 
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Figure 7. Same Figure[6]for the luminosity-dust mass LpiR—Md relation. 
The dashed lines are tracks in constant dust temperature with (left to right) 
T d = 20, 30, 40, 50, and 60 K. 



10 Mq. We find the objects in our sample all fall within a rela- 
tively narrow range of 35 ^ Td ^ 45, though they vary by an order 
of magnitude in Lfir- This is similar to both the LZS and IZS, 
but characteristically different from the SMS, which resembles a 
constant M<j selection with a broad Td distribution extending to 
lower temperatures. The QS, by contrast, does not follow a clear 
trend in Lfir, thought those with high luminosities Z 10 n L© 
tend to have higher dust temperatures and lower dust masses than 
U/HyLIRGs of comparable luminosity. 

This is further illustrated in Figures[7]and[8] where we present 
the Aid — Lfir and Td — Md scalings respectively. In Figure|7] 
the dashed lines represent tracks in constant (left to right) Td = 30, 
40, 50, and 60 K, while in Figure [8] they represent tracks (left to 
right) in constant Lfir = 10 11 , 10 12 , and 10 13 L . These two 
projections also show the dust selection function of the different 
samples samples, whereby the LZS (roughtly), IZS and our sample 
follow a track in constant temperature, the SMS is constant in dust 
mass, and the QS is preferentially warmer and high luminosity with 
lower dust mass. 
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7 THE FAR-IR/RADIO CORRELATION AT HIGH 
REDSHIFT 

It has long been recognized that there exists a tight correlation be- 
tween the far-IR and synchrotron radio emission in s tar-forming 
galax i es, spanning nearly four decades in luminosity ( Helou et al. 
Il985t ICondon & Broderick|[l99ll ; ICondon|[l992l ; lYun et al.ll20oTl) 7 
Though the origin of this correlation is an active subject of on- 
going research, the popular wisdom states that it arises from a 
"cosmic conspiracy" that fixes the ratio of synchrotron losses due 
to magnetic fields associated with massive star formation and 
inverse-Compton from the ambient radiation field dHummelfl 986; 
iHummelet all 19881 ; ICondonll 19921) . This hypothesis is apparently 
born out in observations of local starbursts, which show a roughly 
constant far-IR/radi o correlation over f our decades in magnetic 
field energy density dCondon et alj|l99ll) . 

However, should the far-IR/radio correlation arise from this 
fixed ratio of synchrotron and inverse-Compton losses, one might 
expect it to evolve with redshift. The physics that go vern this ra- 
tio are complex and n ot entirely understood (see e.g., IVoelklll989l ; 
Thom pson et alJ boOfi). Therefore, it is important to check if it re- 
mains valid and/or constant at high redshift and in more luminous 
systems. 

At intermediate redshift 0.3 < z < 1, there is some evi- 
dence for the constacy o f the far-IR/radio corre l ation i Garrettl2002l ; 
iGruppioni etafl 120031; lAppleton et ail |2004 iBovle et all 120071) . 
Howe ver, at high redshift z > 1.5, the picture is somewhat less 
clear jKovacs et ail 120061 : Iviahakis et ail |2007| ; llbar et all 120081 ; 
ISaiina et al.ll2008l) Tand the effects of sample selection - in particu- 
lar, contamination by radio AGN - are uncertain. Our sample offers 
an attractive opportunity to probe the far-IR/radio correlation in a 
relatively narrow redshift range, and for a well-defined sample. 
Overall, our sample has a mean < q >= 2.23 ± 0.04 (see 
f3 — 1.5 and does not vary significa ntly between objects - which is Equation [5} with intrinsic dispersion a q = 0.19 - well in line 

with results for local systems. In fact, only one object (egsb2) is 
inconsistent with the local q, and it is known to harbor an X-ray 
luminous AGN (Huang et al. in prepration). This supports the con- 
clusion of Huang et al. (2008, in preparation) that AGNs contribute 
little to the bolometric luminosity of these objects. Furthermor e, it 
is substantially different from the resutls of Sajina et al. U2008h for 
obscured 24/im sources, which cover a wide range in q both for 



Figure 8. Same Figure[6]for the luminosity-dust mass LpiR—Md relation. 
The dashed lines are tracks in constant far-IR luminosity with (left to right) 
Lfir = 10 11 , 10 12 , and 10 13 L m . 



ity Lfir can be related to the dust mass Md and temperature Td 
via the dust mass absorption coefficient k\ dDe Breuck et al .120031 ; 
I Yang et alj|2007l) : 



47rM d / n\B v (\, Td)dv 



(6) 



where L5 V is the thermal dust spectrum at rest-frame wavelength A. 
For greybody thermal emission, this yields: 



Li 



8uh X^kx fkT d \ 



4+/3 



r(4 + £)C(4 + /3)M d 



(V) 



V h J 

where fi is the dust emissivity, F is the gamma function, and ( is the 
Riemann Zeta function. As a toy model, we use the Milky Way dust 
model of IWeingartner & Draind (2001) at a reference wavelength 
of 125/im (Ki^Sfim ~ 17 g _1 cm 2 ), and make the assumption that 



approximately the case locally (e.g., Dunne et al.|2000l ; lKlaas et"al] 



|2001|; lYun & Carillill2002l ; lYang & Phillip j|2007h . Also, for all the 
following comparisons the total far-IR luminosity is taken from 
40-1000/im (Lfir = £40-1000^111) unless otherwise stated. 

In Figure [6] we show the Lfir — Td relation for our sources 
as compared to the LZS, IZS, SMS, and QS. The dashed lines rep- 
resent tracks (left to right) in constant Md = 10 7 , 10 s , 10 9 and 
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Figure 9. The evolution of the far-IR/radio correlation with redshift param- 
eterized by q (see Equation [5] and §[5]for definitions). The labeling is the 
same as Figure [6] and the grey shaded regio n indicates the local value and 
associated dispersion from lYun et alj (2001). 
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AGN (see Figure 5 of iRisaliti elaLl llOOO). Also, while egs24 is 
not detected in the X-rays, its 25-60^tm flux density ratio suggests 
the presence of a Compton-thick AGN (N H > 4 x 10 24 cm -2 ; 
IRisaliti et alfeOOol) . The remaining 8 objects in our sample have a 
stacked X-ray flux corresponding to < Lx >~ 10 42 ergs s -1 , 
along with < L IR >= 5 x 10 12 L Q (< L x > / < L IR >;£ 
5 x 10~ 4 ) and < L2B^ m / '-^eo^m >~ 0.15, which are all consis- 
tent with a pure starburst. 

In Figures|9]and[T0]we present the q scaling with redshift and 
far-IR luminosity Lfir for our sample ( red circles), loca l systems 
(grey shaded region with q = 2.34 ± 0.2; lYun et al.l200lh . the IZS 
(open stars), and the SMS (black squares). Immediately it appears 
that the local far-IR/radio correlation applies over a wide redshift 
range 0.2 £ z 3 over more than two decades in luminosity. 
The apparent inverse trend of q with both Lfir and redshift in 
the IZS is probably due to their sample selection, which is more 
likely to pick out IR bright AGN at higher redshifts (and luminosi- 
ties). By contrast, both our sample and the SMS show a clear trend 
towards higher q at higher Lfir'- the Spearman rank correlation 
coefficients are p = 0.52 at ~ 2a and p = 0.79 at ~ 3<r for our 
sample and the SMS respectively. This trend is not seen in surveys 
of loc al systems spanning 10 10 Lq ;£ Lfir & 10 12 Lq dYun et all 
1200 lb . However, it is not clear whether or not this effect is due to 
sample selection - our sample and the SMS are selected in very dif- 
ferent ways, and show different trends in dust properties (see § [6]l 
- or the strong redshift evolution of the bright end of the galaxy 
IR luminosity function dSanders & Mi rabel 19961; ILe Floc'h et al.l 
20051: iHopkins et alj|2 008b) and that of massive galaxies (e.g., 
Consel ice et al l2007bb" versus t he comparatively weaker e volution 



of flat-spectrum radio sources jPunlop & Peacock| [l990). Such a 
selection effect would potentially bias z £ 1.5 systems with lumi- 
nosities of Lfir ~ 10 10 ~ n to radio-loud AGN as star formation 
moves preferentially to higher luminosity objects. 
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Figure 10. The correlation between the far-IR/radio correlation coeffi- 
cient q (see Equation \5\ and §|5]for definitions) with the far-IR luminos- 
ity Lfir = ^40— lOOOMm- The labeling is the same as Figure|6] and the 
grey shaded re gion indicates the local value and associated dispersion from 
lYunetaljj200ll) . 

strong and weak PAH emitters - the mean < q >= 1.90 ± 0.13 
with a q = 0.34 for their "strong-PAH" sources with radio de- 
tections, and < q >= 1.38 ± 0.10 with a q = 0.47 for their 
"weak-PAH" sources. This is, however, not surprising consider- 
ing th is sample selection is k nown to be biased towards obscured 
AGN dWeedman et"ail 2006b). We take this as evidence to suggest 
that our selection criteria much more efficiently picks starburst- 
dominated systems in a narrow redshift range. 

The radio/far-IR properties of this sample are also consistent 
with the AGN contribution inferred from their X-ray/far-IR prop- 
erties. The rest-frame 25-60/im flux density ratio versus Lx / Lir 
- where Lx is the rest-frame full band 2-10 keV X-ray luminos- 
ity - has been shown to be sensi tive to AGN activity dde Grijp et al.l 
1 198511 19871 : IRisaliti et al.ll2OO0b . Huang et al. (in preparation) find 
that only three objects in this sample have even marginal 2a X- 
ray detections: egsl (Lx ~ 10 42 ergs s _1 , Lx/Lir ~ 10 -3 , 
L25^ m /L (i0tim « 0.2), egsl4 (Lx ~ 10 42 ergs s~\ L x /Lir as 
3 x lO" 4 , L 2Blim /L 60tim w 0.06), and egsb2 (L x » 10 43 ergs 
s -1 , Lx/Lir tn 10~ 2 , L25 M m/i60 M m ~ 0.15). Of these, only 
egsb2 is consistent with a significant contribution from an obscured 



8 CONCLUSION 

We present far-IR observations, including 1200/xm photometry 
with MAMBO, of high redshift ULIRGs selected by Huang et al. 
(in preparation) on the basis of their observed mid-IR colour. All 
the objects in the sample have mid-IR spectroscopic redshift and 
lie in a remarkably tight redshift window centered at z ~ 2. We 
successfully detected (S/N > 3) 9 of 12 total sources. When com- 
bined with existing MIPS photometry at 70 and 160/im, this pro- 
vides constraints on both sides of the cold dust peak, allowing a 
more robust determination of the far-IR SED. 

We fit a model to the photometry, from which we estimate 
both the far-IR luminosity and dust properties - including tem- 
perature and mass. When compare to other samples, including IR 
luminous systems at low and intermediate redshift, submillimetre 
selected galaxies, and QSOs, and find that our sample, along with 
the low and intermediate redshift sources, is roughly a constant dust 
temperature selection. Submillimetre sources, by contrast, are more 
similar to a constant dust mass selection. Finally, we use existing 
radio observations at 20cm to test the far-IR/radio correlation at 
high redshift. We find that all but one of the sources in our sample 
is consistent with the local relation, suggesting that it remains valid 
at high redshift. This furthermore indicates that the Huang et al. 
(in preparation) method efficiently selects massive, high redshift, 
starburst-dominated systems in a remarkably tight redshift range. 
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